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The magnon energy and amplitude renormalization due to 
intraband particle-hole excitations are studied in a metallic 
antiferromagnet. The change in sign of the intraband contri- 
bution with Lo results in significant differences between static 
and dynamical behaviours. For electron doping, while the 
coherent magnon peak looses spectral weight and is shifted 
to higher energy with doping, the low-energy incoherent part 
becomes increasingly prominent as it narrows in width and 
shifts to lower energy. Implications for spin dynamics in the 
electron doped cuprate Nd2-xCexCu04 are discussed. Due 
to an exact cancellation of two logarithmically divergent spin- 
fluctuation processes, the AF order parameter is unaffected 
to leading order. For hole doping, short wavelength trans- 
verse perturbations (g > g* ~ ^a;) are found to be stable, 
implying short-range AF order with a spin correlation length 
~ l^/x. 
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I. INTRODUCTION 

Antiferromagnetism in the doped cuprates strongly de- 
pends on the type of doping. While the electron doped 
cuprate Nd2-xCexCu04 retains AF order up to a doping 
concentration of about 15%, [|]-|| only 2% hole concen- 
tration destroys AF order in La2_xSrxCu04. 

The simplest microscopic description of this strong 
doping asymmetry emerges within the t — t' Hubbard 
model, with nearest-neighbour (NN) and next-nearest- 
neighbour (NNN) hopping terms t and t', respectively. 
For negative t' (as usually assumed for cuprates) and 
electron doping, the AF state is stable for a range of 
doping concentration, whereas the AF state becomes un- 
stable with respect to transverse perturbations in the AF 
order for any finite hole doping. |^,|| Recently the mag- 
netic phase diagram of the doped t — t' Hubbard model 
has been obtained in the t' — U space, showing the various 
regions of stability and instability with respect to both 
longitudinal and transverse perturbations in the AF or- 
der. D 

The appropriate Hubbard model parameters for 
La2Cu04, determined recently by fitting the spin- wave 
dispersion obtained from high resolution inelastic neu- 
tron scattering studies, |^|-^ indicate that U/t 8. From 
the additional fact that the critical electron doping con- 
centration in Nd2-xCexCu04 is around 20%, the t' — U 
phase diagram then yields t' /t ^ 0.25. This value of t' 
falls in the range 0.15 to 0.5 estimated from band struc- 
ture studies, photoemission data and neutron-scattering 
measurements of high-Tc and related materials. iPHl^ 



A stable AF state of the doped t — t' Hubbard model 
provides a microscopic realization of a metallic antifer- 
romagnet, in which the Fermi energy lies within a quasi- 
particle band. In this paper we will quantitatively study 
the interplay of spin fluctuations and particle-hole excita- 
tions in the metallic AF, focussing on the role of the intra- 
band particle-hole excitations on the spin excitation spec- 
trum, magnon renormalization, and spin-fluctuation cor- 
rection to the AF order parameter. These results should 
be particularly relevant to the electron doped cuprate 
Nd2_xCexCu04, in which metallic conductivity result- 
ing from electron doping suggests that the electrons are 
mobile, jl^ and doping with Ce clearly reduces the spin- 
stiffness constant from the value in the undoped system. 
1^ Metallic antiferromagnetism has also been reported in 
K- (BEDT-TTF)2X, Q V2_x03, @ and NiS2-xSex. 

The metallic antiferromagnet is characterized by 
magnon decay into intraband particle-hole excitations, 
and the resulting magnon damping (F) was studied re- 
cently in the long wavelength limit (g << 1). |^ In terms 
of the Fermi circle radius a = y/2TTX of the doping pockets 
formed around (±7r, 0) and (0, ±7r), where x is the doping 
concentration, a new doping-dependent energy scale At'a 
was identified, the relative magnitude of which, in com- 
parison with the magnon energy scale \/2J, should essen- 
tially determine the spin excitation spectrum. The imag- 
inary part of x°(Q5 ^) was found to vanish for uj > At'aq, 
and increase linearly with to for to << At'aq. This im- 
plies that the magnon spectrum should exhibit a sharp 
peak at « \f2Jq for \f2J > 4i'a (low doping limit), 
and magnon broadening should appear with increasing 
doping when \/2J < At' a. 

In addition to the spin excitation spectrum for all q 
(section HI), we will also study the quantum correction 
to sublattice magnetization (section IV), to determine 
whether long-range AF order survives in the metallic 
antiferromagnet when quantum spin fluctuations are in- 
cluded. At the one-loop level, the spin-flip process is ac- 
companied by virtual magnon emission and absorption, 
and therefore we will also examine whether the magnon 
propagator (whose amplitude goes like l/g in the AF in- 
sulator) is significantly renormalized due to the intraband 
particle- hole excitations (section II). The consequences of 
magnon renormalization on finite temperature spin dy- 
namics and Neel temperature in electron doped cuprates 
are briefiy discussed in section V. 

We consider the t — t' Hubbard model on a square 
lattice, with NN and NNN hopping terms t and t' con- 
necting sites i to i + S and i + n, respectively: 
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NN 



NNN 



(1) 

In the following we set t — 1. Since a particle-hole trans- 
formation maps the t' model with hole (electron) doping 
on the —t' model with electron (hole) doping, the pos- 
itive (negative) t' model and hole (electron) doping is 
appropriate to study for the electron-doped compound 
Nd2-xCexCu04, for which a negative t' is usually as- 
sumed. 



II. MAGNON RENORMALIZATION 

In this section we consider positive t' and hole doping, 
so that the results are appropriate for the electron-doped 
cuprates. We also consider the strong coupling limit for 
analytical simplicity. To study transverse spin fluctua- 
tions in the doped antiferromagnetic state we evaluate 
the time-ordered magnon propagator 

X-+(q,^)= / dt5]e'(-*-'i-"---)(vI/G|T[^-(t)5+(0)]|*G) 

(2) 

involving the spin-lowering and spin-raising operators 
and Sj' at lattice sites i and j. In the random phase 
approximation (RPA), we have 



x"(q,^) 

1-Ux°{ci,u;)- 



(3) 



where x'^(q, w) is the zeroth-order, antiparallel-spin 
particle- hole propagator, evaluated in the broken- 
symmetry, Hartree-Fock (HF) state of the metallic anti- 
ferromagnet, and involves both intcrband and intraband 
excitations. We examine the q, uj dependence of x'^ili 
which essentially determines the magnon energy and am- 
plitude renormalization. We focus on the real part of 
[x^lq^^)] ill the underdoped limit, as its imaginary part 
was shown to vanish for u > At'aq in the long wavelength 
limit. Hi 

The HF-level description of the metallic AF state 
has been discussed earlier. |^ The NNN hopping term 
modifies only the AF-state quasiparticle energies i?k, 
but not the quasiparticle amplitudes. We have = 
e((. ± + e^, for the upper and lower Hubbard 

bands, where ek = —2t{coskx + cosky) and ej^. = 
—At' cos kx cos ky are the free-fermion energies corre- 
sponding to the NN and NNN hopping terms. Here 
2A = mil in terms of the sublattice magnetization m. 

In the two-sublattice basis (labelled by A and B), the 
intraband contributions to the diagonal and off-diagonal 
matrix elements of [x''(q, w)] are obtained as: 
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= [x°(q,^)]l?r 
(4) 



where X^k indicates that states k are below the Fermi 
energy Ep, while states k — q are above Ep, and Aq = 
E^-q — -Bk is the intraband particle-hole energy differ- 
ence in the lower Hubbard band. 

The above three intraband contributions can be writ- 
ten, in units of (i^/A^), as {aq — bqUj/2J), (aq + 6qti^/2J), 
and Cq, respectively. The dimensionless coefficients 
flqj^qiCq arc all functions of cj^ and, for u ^ 0, vanish 
quadratically with g as g — > due to phase space re- 
striction. However, for oj oc q (near the magnon-mode 
energy), they all have well defined limits as g — > 0. 

Including the interband contributions as well, which 
were studied earlier for finite doping up to order x^, j|] 
[x'^(q, a;)] can be written as 
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The term represents modifications in the interband 

contribution due to NNN exchange energy J' — At jU 
and finite doping, and to first order in a;, is given by 
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(6) 



where 7q = (cos^a, + con qy)/ 2 and 7^ — cos q^ cos qy). 
The scaled frequency w = m'^ {l+bq){uj /2J) = /q (w/2J), 
where m = 1 — x, reflects the frequency-scale renormal- 
ization due to doping and intraband excitations. 

Substituting x''(q, w) from Eq. (5) in Eq. (3), the 
RPA-level magnon propagator is obtained as 



(1 - "q - aq) 
, -(7q-Cq) 

2J 
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where 



r!q(c^) = 2J/qi[(l 



(7q 



x2ll/2 



(7) 



(8) 
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is frequency dependent due to the intraband contribu- 
tion, resulting in the magnon amphtude renormahzation. 
The magnon-mode energies for momentum q are given by 
the poles at ujq = ±f2q(wq). 

For J << 8t', as in the strong coupling limit, the 
magnon renormalization due to bq is negligible, and in 
the following we take the frequency-scale renormalization 
factor /q — m^(l + bq) w 1, for simplicity. 

In the long- wavelength and low-doping limits (q << 
1, X << 1) we have 7q 1 - q'^/A, 7q « 1 - and 
J, Cq << 1. Substituting for a' from Eq. (6) in Eq. 



(8), the expression for Qq(uj) simplifies to 

ftqiuj) = V2J[aintcr - ftintra (w)] "^^^^ 

where the dimensionless quantities 



Q^intcr 
and aintra{L0) 



1- {2J'/J){1- Ax/t') 

4(aq - Cq)/q^ 



(9) 



(10) 
(11) 



were introduced earlier as coefhcients of the q^ term in 
the interband and intraband contributions to the trans- 
verse response eigenvalue in the context of the stability 
analysis of the doped AF state of the t ~ t' Hubbard 
model. 1^ From Eq. (4) which defines Cq, 6q, and Cq, 
and substituting for ek-q — Ck « 2tqa cos 0, and the intra- 
band particle-hole energy difference £^k-q " = ~ 
At' qa cos 9 in the limits t' > Jx and g << 1, we obtain 



Ax 

aintra(,W) = — 



cos^ - {oj/ At' qay 



(12) 



As {uj) is a function of uj/q, it follows from Eq. (9) 

that the magnon-mode energy ujq is proportional to q, as 
in the AF insulator. 

The change in sign of the intraband coefRcient 
aintra('^) at uj/At'qa = 1 (see Fig. 1) has interesting con- 
sequences on the doping dependence of the magnon-mode 
energy. In the static hmit, aintra(O) — Ax/2t' is posi- 
tive, and the intraband contribution softens the magnon 
mode, eventually leading to the instability of the AF 
state when aintra = Winter • @ However, for lu > At'qa, the 
intraband coefficient aintra is negative, implying that the 
intraband excitation actually stiffens the magnon-mode. 

Therefore, for V^J > At'a = At' \J2'kx (underdoped 
limit) , the coherent magnon peak at LOq w ^/2Jq is shifted 
to higher energy, this shift being proportional to x^ in the 
small doping limit. With increasing doping, when a/2 J < 
4i'a, the magnon peak is broadened due to damping and 
gets shifted to lower energy due to the change in sign of 
Q;intra('^)- As tlic magnou damping term goes as u for 
to « At'qa, the magnon peak significantly narrows 
down as it shifts to lower energies with doping. 

Doping stiffens the magnon mode through the inter- 
band contribution as well, by suppressing the softening 
due to the J'-induced frustration. This is clearly seen 
from the expression for the interband coefRcient aintor in 
Eq. (10), which increases with x. 




FIG. 1. The intraband coefficient Qfintra('jj) changes sign 
&t u) — At'qa, implying that the magnon-mode softening 
for ujq < At'qa changes into magnon-mode stiffening for 
i^q > At'qa. 

Turning now to the magnon amplitude renormal- 
ization, by expanding the frequency-dependent energy 
flq{uj) near the mode energy cjq, we obtain 



to — ilq(a') 



duj 



{uj - LUq) 



(13) 



which yields the magnon renormalization factor 

, -1 



Z 



1 
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(14) 



From Eqs. (9) and (12), we obtain 



V2j\ Ax ^0 



At'a 



cos 6 



t' 7_W2 TT [cos^ 9- {LOq/ At'qa)^]^ 



(15) 



which is negative, and vanishes with the doping concen- 
tration as x^, through the x dependence of the integral. 
These consequences of doping, namely the stiffening of 
the magnon mode and reduction in the coherent spectral 
weight for low doping, changing into magnon softening 
and broadening with increasing doping, are confirmed in 
a numerical study of the magnon spectral function, as 
discussed below. 



III. SPIN FLUCTUATION SPECTRUM 

The spectrum of transverse spin fluctuations in the 
doped AF state, evaluated numerically from the imag- 
inary part of the magnon propagator x (q, u) , is dis- 
cussed in this section. In terms of the two complex eigen- 
values A„(q, w) of the x°(q, lj) matrix, we obtain for the 
spin fluctuation spectral function Aq{uj) 
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FIG. 2. The magnon spectral func- 

tion Aq{uj) = Tr Im x '"(li'^) for long wavelength mode 
q — (0.2,0), showing the shift of the coherent magnon peak 
to higher energy with doping, accompanied by a reduction of 
the magnon amplitude. With increasing doping the broad in- 
coherent spectrum acquires strength and narrows into a peak 
which shifts down in energy. 




FIG. 3. For the short-wavelength mode q = (1,0), the 
magnon peak broadens almost immediately with doping, with 
no change in the peak energy. With increasing doping, the 
magnon peak becomes strongly asymmetric, rapidly narrow- 
ing into a strong peak at low energy. 



A,i(q,w) 



n=l,2 



- UXn{q,Uj) 



(16) 



In the numerical evaluation of the x'^(q, w) matrix, we 
have taken a grid size dk^ = dky = 0.01 and an imaginary 
term 77 = 0.001 throughout. The spectral function Aq{uj) 
has been recently studied for the AF insulator in the full 
U range from weak coupling to strong coupling, includ- 
ing contributions from both single-particle and collective 
(magnon) excitations. |17| In the following, we discuss 
the results in the context of cuprates, separately consid- 
ering electron doping (positive t' and hole doping) and 
hole doping (negative t' and hole doping). 



FIG. 4. For a fixed doping concentration, the magnon 
spectral function for different q shows a linear dependence of 
the coherent magnon peak energy on q. With decreasing q 
the incoherent part is compressed into a narrow peak, which 
competes in strength with the coherent peak. 



A. Electron doping 

The results for Aq{uj) are shown in Figs. 2-5, 
which confirm the characteristic features obtained in the 
magnon renormalization study of section II in the small 
q limit. With doping, the coherent magnon peak shifts 
to higher energy with decreasing amplitude (see Fig. 
1). Simultaneously, the incoherent part of the magnon 
spectral function at low frequency progressively becomes 
narrower, shifts to lower energy and develops oscillator 
strength. The softening is due to the intraband contri- 
bution and the narrowing follows from the behaviour of 
Im which decreases with lu. 

On the other hand, for short wavelength modes, the 
magnon peak is seen to broaden almost immediately with 
doping (see Fig. 3), with no change in the peak energy. 
A study of the imaginary part of x°(q, w) shows that 
the energy range over which Imx'^(q, (jj) is finite initially 
increases as 4,1' aq, but then increases more rapidly with q, 
eventually overtaking the magnon energy. Therefore for 
short wavelength modes, the magnon energy falls within 
this range, resulting in the broadening. 

The magnon softening is seen to be moderately q de- 
pendent. While the magnon peak shifts to nearly a; « 
at a; ~ 10% for small g, a substantial energy gap still 
remains for g ~ 1 (see Fig. 3). This interestingly shows 
that the instability of the AF state with respect to trans- 
verse fluctuations in the order parameter strongly de- 
pends on the fluctuation wavelength. While the AF state 
becomes unstable with respect to long wavelength modes 
aX X = Xc, the short wavelength modes are still positive 
energy modes, so that short-range AF order should sur- 
vive. 



4 





1 1 1 


1 1 


ii 9. = 0.5-^ 
















!i 1.0 _ 








!i 1.5 - 


_ 






!i 2.0 

1 j 








j j 

i U ^8 








\ \\ *' = 0.2 - 








\. 1 j x = 7% - 








Mi q=(<?.,0)- 












1 . .1,. 


t ^ 1 


■ 1 . ■ 1 1 



0.1 0.2 0.3 0.4 0..5 0.6 0.7 0.8 0.9 1 




FIG. 5. The incoherent part of the magnon spectral 
function decreases progressively with increasing wavevector 
q — (<7rc,0), for a fixed doping concentration of a; = 7%. 
Magnon damping and linewidth sharply decrease as ap- 
proaches ~ 2.0, and for > 2 only the coherent magnon 
peak remains at low energy, with an insignificant incoherent 
part at higher energy (not shown). 

For a fixed doping concentration, the coherent magnon 
peak energy is proportional to the wavevector q in the 
small q limit (see Fig. 4). For short wavelength modes, 
the incoherent part decreases dramatically with increas- 
ing wavevector (see Fig. 5). Magnon damping and 
linewidth sharply decrease as q approaches ^ 2.0, and 
for q > 2 only the coherent magnon peak remains at low 
energy, with an insignificant incoherent part at higher 
energy. 

As the doping approaches the critical concentration Xc, 
above which the AF state is unstable, the appearance of 
narrow magnon modes at very low energy despite their 
strong renormalization due to the intraband particle-hole 
excitations, is a noteworthy feature of the metallic AF 
state. 



B. Hole doping 

The strong intraband contribution renders the homo- 
geneous AF state unstable for any amount of hole dop- 
ing. This instability with respect to transverse per- 
turbations in the AF order is signalled by the trans- 
verse response eigenvalue UXq exceeding unity for long 
wavelength (small q) modes, indicating absence of long- 
range AF order, and signalling a tendency towards in- 
commensurate ordering with wavevector different from 
Q = (7r,7r). Different types of homogeneous spiral 
phases, and their stability with respect to longitudinal 
and transverse perturbations have been studied in detail. 

i 

However, short-range AF order appears stable, as in- 
dicated by the full q dependence of UXq (see Fig. 6) for 
the two eigenvalues Aq of the x°(q) matrix for q — {q, q) 
in the range < g < tt. The reflection symmetry about 



FIG. 6. The q dependence of L'^Aq in the range < g < tt 
for the two eigenvalues Aq of the x"(ci) rnatrix for hole con- 
centration of 5%. 
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FIG. 7. The g dependence of UXq for different hole doping 
concentrations x, showing that UXq drops below 1 above some 
value q = q* . Note that the small-q parabolic behaviour of 
UXq is independent of doping concentration. 

q = 7r/2 is due to the location of the hole pockets around 
(±7r/2, ±7r/2) in the Brillouin zone. Considering the 
larger eigenvalue for small q, the parabolic behaviour of 
UXq is independent of doping concentration (see Fig. 7), 
as obtained earlier. However, beyond a certain q value 
U Aq starts decreasing and eventually drops below unity 
at g = q* . This signals the stability of the AF state 
with respect to transverse perturbations of wavelength 
shorter than A* = 2TT/q*, implying that short-range AF 
order can exist in the hole doped AF upto length scale 
^ ^ A*. For X > 15%, the transverse response eigenvalue 
UXq > 1 for all q, indicating no stability for AF domains 
of any size. 

The q dependence of UXq for different hole doping con- 
centrations (Fig. 7) shows that the wavevector value 
q* where UXq drops below unity increases with doping 
concentration x, indicating diminishing spin correlation 
length. The behaviour of the spin correlation length 
^/a = 2TT/q* with x is shown in Fig. 8, clearly show- 
ing a ^l\fx dependence, as seen in neutron scattering 
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experiments. 
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FIG. 8. Plot of {q*)^ vs. X shows a linear dependence, 
implying a x~^^'^ dependence of the spin correlation length ^ 
on hole doping concentration x. 
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FIG. 9. The imaginary part of the particle-hole propagator 
■X^{<\,ijj) shows a linear behaviour over a large uj range. 

Figure 9 shows that imaginary part Im[x°(q, w)]aa on 
the A sublattice has a remarkably linear behaviour over 
a large u range. A similar behaviour for q not too small 
is also observed in the electron-doped case. This low- 
frequency contribution in Imx"(q, to), arising from the in- 
traband excitations in the antiparallel-spin particle-hole 
propagator x"(q, ^^), is suppressed to order /U"^ due to 
the AF coherence factors. However, no such suppression 
is present in the parallel-spin particle-hole propagators 
[7r°(q, Ci')]AA or [7r°(q, w)]bb, where both particle and hole 
amplitudes (in the lower band) are of order unity. Such 
parallel-spin particle-hole propagators are involved in the 
charge density fluctuations. 

Although the AF state is unstable at the static level, 
the spectral function shows sharp peaks at finite fre- 
quency (see fig. 10). This is due to the change in 
sign of the intraband coefficient aintra(w) with increas- 
ing oj, discernible from the change in curvature of the 
spectral function. As for the electron-doped case, the 
small-g magnon peaks are sharp because Imx"(q, ljj) van- 



ishes at the magnon-mode energy, whereas at higher q the 
magnon-mode energy lies within the energy scale over 
which Imx° (q, w) is finite 




0.50 



FIG. 10. The spectral function shows sharp magnon peaks 
for low q, and broadened peaks at higher q due to magnon 
damping. 
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FIG. 11. For fluctuation modes of wavelength smaller than 
the domain size, {q > q* ~ 0.9 for x = 0.05), the magnon 
modes are strongly softened, with a small incoherent part 
appearing at higher energy. 

(see Fig. 9), resulting in magnon broadening. 

Figure 11 shows that there is an abrupt change in the 
spectral function when q crosses q* w 0.9 for x — 0.05 
(see Figure 7). For fluctuation modes of wavelength 
smaller than the AF domain size {q > q*), the magnon 
modes are strongly softened by the intraband contribu- 
tion, with a small incoherent part appearing at higher en- 
ergy. In contrast, in the electron-doped case, the coherent 
part was shifted to higher energy and a weak incoherent 
part appeared at low energy. Both magnon damping and 
softening have been observed in the hole-doped cuprate 
La2_xBaxCu04. [fl 
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IV. FINITE TEMPERATURE SPIN DYNAMICS 

At finite temperature T << J, the decrease in sub- 
lattice magnetization m{T) due to thermal excitation 
of long wavelength magnons has been studied in a 
highly anisotropic layered antiferromagnet with planar 
exchange 
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FIG. 12. The long-wavelength magnon energy scale 
^poak^^ falls linearly with the doping concentration x, imply- 
ing a corresponding decrease in the spin-wave stiffness con- 
stant and the Neel temperature. 

interaction energy J = At'^ /U, an effective interlayer hop- 
ping tz , and an effective anisotropy gap Adm due to 
the Dzyaloshinski-Moriya interaction. An estimate of 
the Neel temperature Tn was obtained from the predom- 
inant T In T faU off of m(T) with T as fcBlN ~ J/ ln(l/r) 
for Adm << 2Jr, and /cbTn ^ J/ln(J/ADM) for 
2Jr << Admj where r — tz/t is the interlayer-to-planar 
hopping ratio. In both cases, Tn is proportional to the 
long wavelength magnon energy scale (wq/g)g_>o ^ J- 

Figure 2 shows that the relevant magnon energy scale 
which will essentially determine the low temperature spin 
dynamics in the doped AF is the incoherent magnon peak 
energy, which shifts to lower energy and gains spectral 
weight with increasing doping. The doping dependence 
of this magnon energy scale uj^°^^/q shows a nearly linear 
decrease with doping concentration (see Fig. 12). 

This nearly linear decrease in the incoherent magnon 
peak energy with doping concentration provides an ex- 
planation for the observed decrease in unison of the 
spin-wave stiffness constant ps and the Neel tempera- 
ture Tn with Ce doping in both Nd2-xCexCu04 and 
Pr2-xCexCu04. A spin dilution model, within which 
a Cu^"*" S=l/2 ion is assumed to be converted into a lo- 
calized Cu^+ S=0 ion for each added Ce ion, has been 
found to quantitatively describe the x-dependence of ps 
and Tn- However, the localization of the added elec- 
trons in this scenario contradicts the metallic conductiv- 
ity resulting from Ce doping. 



V. SPIN-FLUCTUATION CORRECTION TO 
SUBLATTICE MAGNETIZATION 

Figure 13 shows the first-order spin-fluctuation correc- 
tion to the electron propagator. Physically, the process 1 
(2) represents a spin-t electron in state k below (above) 
the Fermi energy emitting a virtual magnon with mo- 
mentum Q and flipping into a spin-| electron in state 
Q,0 



k, oj 
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t 



k 



Q 

■ 



k, oj 
t 



FIG. 13. The quantum spin-fluctuation correction to the 
electron propagator, representing a spin flip accompanied by 
the emission and absorption of a virtual magnon. 

k — Q above (below) the Fermi energy, which then re- 
absorbs the virtual magnon. Such spin-flip processes in- 
volve transfer of spectral weight across the Fermi energy, 
and therefore result in corrections to electron densities 
7T,| and ni, and hence to the sublattice magnetization. 

Considering first the process (1) involving < E-p 
and £^k-Q > Ep, the integrated spectral weight trans- 
ferred above the Fermi energy yields the reduction in the 
spin-l density on an A-sublattice site 



-6'i 



.(1) 



uj)\aa 



= ^ r^Im[5GT(k, 
1, J ^ 



a 



k~Q,i 



{Ei,-Q -Eu + n)^ 

(17) 



involving the retarded part of the magnon propagator 
XR^(Q,fi)- Here Ok-Q i ^he spin-| electron density 
on the A sublattice site, and in the strong coupling limit, 
is given by 1 and e^_q/4A^ for state k — Q in the upper 
and lower band, respectively. 

For the undoped antiferromagnet, where the spin-fiip 
process involves an interband excitation, and states k and 
k — Q are energetically separated by the AF gap 2 A, this 
correction reduces the A-sublattice density from 1 to 
0.8 for the square lattice. |2^ A similar process transfers 
spin- J, spectral weight from the upper (unoccupied) band 
to the lower band, increasing from to 0.2, resulting 
in a net reduction of 0.4 in the sublattice magnetization. 

In the metallic antiferromagnet, however, spin-flip pro- 
cesses involve the intraband excitations as well, and the 
small energy denominator E^ ^ — E^ ~ At'Qa in the 
Q << 1 limit drastically increases the spin-fluctuation 
contribution to the particle density correction. On the 
other hand, the phase-space restriction on states k and 
k — Q to lie across the Fermi energy suppresses the in- 
traband contribution in the long-wavelength limit by a 
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factor Q. In the following, we examine the combined ef- 
fect of these two features of intraband excitations on the 
quantum correction to sublattice magnetization. 

As discussed earlier, except for a renormalized magnon 
velocity c and amplitude Z, the coherent part of the 
magnon propagator for Q « 1 is qualitatively un- 
changed in the weak doping limit {-/2J > 4i'a), with 
Im X]^^(Q,w) ^ {Z/Q)5{ijj — cQ) as for the AF insula- 
tor. Substituting in Eq. (17), we obtain 



5n 



(1) 



dn z 



27r 

5{VL - cQ) 



n/2 
-1x12 



Q dQ r'^ ad0Q cose , 
1 



{At'aQ cosd + n)^ 



(18) 



After doing the il integral, the resulting Q integral in 
Eq. (18) is of the form J dQ/Q, which yields a logarith- 
mically divergent reduction in the A-sublattice density 
n| from process (1). Similarly, the process (2) involv- 
ing > Ep and E^_^ < Ep, and the advanced part 

of the magnon propagator Xa^{Qj transfers spectral 
weight below the Fermi energy, yielding an identical log- 
arithmically divergent enhancement to Remarkably, 
it is this exact cancellation of two logarithmically diver- 
gent contributions which is responsible for the survival 
of long-range AF order in the weakly (electron) doped 
cuprate. 



VI. CONCLUSIONS 

The study of spin fluctuations in the t ~ t' Hub- 
bard model shows that when doped electrons (holes) are 
added to the bottom (top) of the upper (lower) Hub- 
bard band, the consequent intraband particle-hole exci- 
tations strongly renormalize the magnon propagator in 
the metallic AF state. The change in sign of the intra- 
band coefhcient aintia('^) with a; is a key result, leading 
to differences between static and dynamical behaviours. 
Several features of the magnetic properties of the electron 
and hole doped cuprates are understandable within this 
microscopic model. These include: i) finite (nearly zero) 
critical doping concentration Xc above which long-range 
AF order is destroyed in electron (hole) doped cuprates, 
ii) a linear decrease in the spin-wave stiffness constant 
and Neel temperature with doping concentration, iii) 
magnon broadening and softening in hole doped cuprates, 
iv) enhanced correlations at the dynamical level. 

However, important features of the hole-doped 
cuprates require further investigation. While the t — t' 
Hubbard model (with negative t') indicates instability of 
the AF state and tendency towards incommensurate or- 
dering for any finite hole doping, spin fluctuations in the 
hole-doped cuprates are commensurate for low doping 
and a commensurate-incommensurate transition occurs 
at X « 0.05. The exact mechanism responsible for the 



loss of long-range AF order for x < 0.05 therefore re- 
mains unclear. Interestingly, it is in this doping regime 
that a magnon damping term F ~ T in La2-xSrxCu04 
and La2-xBaxCu04 has been observed at finite temper- 
ature, which accounts for the anomalous nuclear 
spin relaxation rate, resistivity etc. 



16 
17 



Electronic address: avinas@iitk.ac. in 
T. R. Thurston, M. Matsuda, K. Kakurai, K. Yamada, 
Y. Endoh, R. J. Birgeneau, P. M. Gehring, Y. Hidaka, 
M. A. Kastner, T. Murakami, and G. Shirane, Phys. Rev. 
Lett. 65, 263 (1990). 

C. Almason and M. B. Maple, in Chemistry of High- 
Temperature Superconductors, edited by C. N. R. Rao, 
World Scientific, Singapore (1991). 

M. Matsuda, Y. Endoh, K. Yamada, H. Kojima, I. 
Tanaka, R. J. Birgeneau, M. A. Kastner, and G. Shirane, 
Phys. rev. 45, 12548 (1992). 

A. V. Chubukov and K. A. Musaelian, J. Phys.: Condens. 
Matter, 7, 133, (1995). 

A. Singh and H. Ghosh, Phys. Rev. B 65, 134414 (2002). 
R. Coldea, S. M. Hayden, G. Aeppli, T. G. Perring, C. 

D. Frost, T. E. Mason, S.-W.Cheong, and Z. Fisk, Phys. 
Rev. Lett. 86, 5377 (2001). 

N. M. R. Peres and M. A. N. Araujo, Phys. Rev. B 65, 
132404 (2002). 



A. Singh and P. Goswami, Report No. cond-mat/0205117 
(2002). 

P. Benard, L. Chen, and A. -M. S. Tremblay, Phys. Rev. 
B 47, 589 (1993); 

A. VeiUeux, A. Dare, L. Chen, Y. M. Vilk, and A. -M. S. 

Tremblay, Phys. Rev. B 52, 16255 (1995). 

T. Tohyama and S. Maekawa, Phys. Rev. B 49, 3596 

(1993). 

G. Stemmann, C. Pepin, and M. Lavagna, Phys. Rev. B 
50, 4075 (1994). 

O. K. Andersen, A. 1. Liechtenstein, O. Jepsen, and F. 
Paulsen, J. Phys. Chem. Solids 56, 1573 (1995). 
Z. Z. Wang, T. R. Chien, N. P. Ong, J. M. Tarascon, and 
E. Wang, Phys. Rev. B 43, 3020 (1991); S. J. Hagen, 
J. L. Peng, Z. Y. Li, and R. L. Greene, ibid. 43, 13606 
(1991); S. Kohiki, J. Kawai, T. Kamada, S. Hayashi, H. 
Adachi, K. Setsune, and K. Wasa, Physica C 166, 437 
(1990). 

T. Sasaki et al, Solid. St. Commun. 75, 93 (1990). 

S. Carter et al, Phys. Rev. Lett. 67, 3440 (1991); 

W. Bao et al, ibid. 71, 766 (1993). 

S. Sudo, J. Mag. Mag Mat. 114, 57 (199 2). 

A. Singh, Report No. |cond-mat/0006079| (2000). 

T. R. Thurston, R. J. Birgeneau, M. A. Kastner, N. W. 

Preyer, G. Shirane, Y. Fujii, K. Yamada, Y. Endoh, Y. 

Hidaka, and T. Murakami, Phys. Rev. B 40, 4585 (1989). 

G. AeppU, S. M. Hayden, H. A. Mook, Z. Fisk, S-W. 

Cheong, D. Rytz, J. P. Remeika, G. P. Espinosa, and A. 

S. Cooper Phys. Rev. Lett. 62, 2052 (1989). 



8 



A. Singh, Z. Tesanovic, H. Tang, G. Xiao, C. L. Chien, 
and J. C. Wallier, Phys. Rev. Lett. 64, 2571 (1990). 
A. Singh and Z. Tesanovic, Phys. Rev. B 43, 11445 
(1991). 

A. Singh and Z. Tesanovic, Phys. Rev. B 41, 614 (1990); 
ibid. 41, 11 457 (1990). 

S. M. Hayden, G. AeppU, R. Osborn, A. D. Taylor, T. 
G. Perring, S.-W. Cheong, and Z. Fisk, Phys. Rev. Lett. 
67, 3622 (1991). 

B. Keimer, R. J. Birgeneau, A. Cassanho, Y. Endoh, R. 
W. Erwin, M. A. Kastner, and G. Shirane, Phys. Rev. 
Lett. 67, 1930 (1991). 



